Vortex bending and tightly packed vortex lattices in Bose-Einstein condensates 
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We study in detail the structure of the ground state of elongated rotating Bose-Einstein condensates. 
Such ground state is composed of one or more vortex lines which bend even in completely symmetric 
setups. This symmetry breaking allows the condensate to smoothly adapt to rotation and to produce 
tightly packed arrays of vortex lines. The dependence of vortex bending with respect to relevant 
parameters is studied. 
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I. INTRODUCTION 

A stirred coffee climbing up the cup walls, a whirl of 
water spiraling down a drain, or a tornado walking over 
the earth, are natural phenomena involving the rotation 
of a fluid, let it be coffee, water or dusty air. In all of 
them the fluid is subject to a balance between an exter- 
nal force and the centrifugal force, and this balance is 
responsible for a depression of the density along the axis 
of rotation. These structures have no trivial symmetries 
as is evident in the case of tornadoes, which have con- 
tinuously changing bent shapes. None of these structures 
have long lives because viscosity, imperfections and other 
dissipative mechanisms play such an important role that 
rotation cannot be self sustained. 

This behavior contrasts with that of the so called su- 
perfluids, which represent a state of matter of negligible 
viscosity or dry fluid The lack of viscosity allows a 
superfluid to host rotation for long times with little or 
no external intervention. Besides, when a superfluid ro- 
tates it follows a special type of flow which is irrotational, 
V x v = 0, except for a finite number of extended sin- 
gularities, which are called vortex lines. These vortex 
lines, first predicted for 4 He condensates, are the super- 
fluid equivalent of robust whirls and tornadoes. However, 
due to centrifugal forces, the density of the superfluid ||] 
becomes zero along the vortex line something which is 
difficult to achieve in ordinary liquids. 

Until recently only two examples of true superfluids 
were known: the superfluid phases of 4 He and 3 He. In 
both cases, the strength of the many-body interactions 
obscure the signatures of superfluidity and make their 
study difficult. The search for superfluid-type, weakly 
interacting systems has led to the experimental realiza- 
tions of Bose-Einstein condensates (BEC) using dilute 
gases of alkali atoms ||. These condensates are ruled at 
low temperature by a mesoscopic quantum wavefunction, 
■0(r) = \J ' p{v)e ie ( r \ whose phase determines a velocity 
field, v = hW8(r)/m, that may host vortices along the 
lines where 9 is not well defined. The existence of vor- 
tices in BEC has been directly confirmed in two sets of 
experiments. The first one is based on the manipulation 
of the internal degrees of freedom of an almost spherical 
condensate M. In the second one, vortices are produced 



by "mechanical means" in rotating elongated traps with 
transverse asymmetries a geometry that leads to 

intriguing effects on the vortex nucleation process. 
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FIG. 1. Ground state of a condensate in an asymmetric 
trap (e = 0.03) as the one from current experiments for an 
angular speed of Q = 0.55u±. Shownare the vortex cores and 
several transverse sections of the condensate. The gray scale 
indicates the modulus of the wavefunction. We also indicate 
the axis of rotation of the trap. 

In Ref. [Q we find several commonly accepted predic- 
tions about vortices in rotating traps. Namely, (i) a ro- 
tating trap leads to the production of one or more vortices 
of unit topological charge m, being m = (2tt)^ 1 § V6. (ii) 
Vortices appear in finite numbers, with straight shapes, 
forming lattices with p-fold symmetries, (iii) Certain 
critical rotation speeds f2 p < fl p +i should be surpassed 
before each new vortex is nucleated, (iv) Asymmetric 
states without rotational symmetry of any kind are found 
to be energetically unstable, which means that there are 
perturbations with no energetical cost that can destroy 
such states. However, Ref. jjj is based on a variational 
ansatz for the weak coupling limit, which lacks longitudi- 
nal degrees of freedom and implicitly induces the p-fold 
symmetry of the vortex lattice. Recent works abandoned 
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these constraints but either focused on infinitely long con- 
densates or studied almost spherically symmetric traps 

In this paper we study an elongated three-dimensional 
Bose-Einstein condensate subject to rotation. Our main 
result is that stationary vortex lines are bent even in com- 
pletely symmetric setups and that the bending depends 
on the type of trap on the number of atoms of the ex- 
periment. The fact that vortex lines may bend was men- 
tioned in Ref. || where longitudinal excitation modes of 
a vortex line induce vortex bending. However the physics 
in that paper differs drastically from our work, in which 
we consider the ground state and the generation of stable 
and stationary structures. Second, in a complementary 
paper ||lo|| we have presented some evidences in favor of 
that bending and how it could be related to the exper- 
imental results of the ENS ||. In this work we extend 
that work to the analysis of arrays of bent vortices and 
its properties as well as the dependence of the bending 
of the vortices on the geometry of the trap. 

Our work consists on two parts. In Section || we pro- 
vide the mathematical foundations for the search of the 
rotating condensate most favorable configurations. To 
do so we first introduce the mean-field model used to 
describe the condensate together with some basic defini- 
tions and a new energy functional, whose minima are the 
possible ground states of the condensate. In this section 
we also develop an efficient minimization method for this 
functional. 



All these tools are then applied in Sect. Ill to different 
configurations. We start with setups from current exper- 
iments g and show that in elongated traps with many 
bosons the vortex lines are bent and form Abrikosov lat- 
tices which are regular only in the core of the gas cloud 
and deform close to the boundaries. Next we study the 
dependence of the vortex bending with respect to the free 
parameters in our theory. Finally in Sect. ^ we offer 
our conclusions together with some open questions. 



II. MATHEMATICAL TOOLS 

A. The mean-field model 

For current experiments it is an accurate approxima- 
tion to use the zero temperature mean field theory of the 
condensate, in which the atomic cloud is described by 
a single wavefunction ip(r,t) ruled by a Gross-Pitaevskii 
equation (GPE). In the case of rotating systems it is use- 
ful to consider the problem on the mobile reference frame 
that moves with the trap |Tl] ] , in which the equation reads 



~A + V (r)+g\ib\ 2 -fLL z 



(1) 



Here L z = i {xd y — yd x ) is the hermitian operator that 
represents the angular momentum along the z-axis and 
the effective trapping potential in given by 



Vb(r) = -Ji 



(1 - e)x 2 + (1 + e)y 2 ] + -uj 2 z z 2 



(2) 



In Eq. (Q) we have applied a convenient adimensional- 
ization which uses the harmonic oscillator length, a±_ — 
yJh/mRbU)^ and period, r = u;J . With these units the 
nonlinear parameter becomes g = Airas/a±. 



B. Variational formulation 

There are several conserved quantities associated to 
Eq. (0). The first one is the norm of the wavefunc- 
tion, N[tp] = J dr, which is related to the number of 
bosons in the condensate. The second conserved equation 
is the energy of the gaseous condensate 



E[fl>] = J ^ 

= E Q {il } )-QL z {ip) 



l -A + V Q {v) + 9 -\iP\ 2 -m z 



ipdr 



(3) 



In current experiments with stirred condensates 
P,p|p^] , the gaseous cloud reaches certain long lived con- 
figurations with one ore more vortices. In this work we 
are not interested on the precise dynamics of the nucle- 
ation process, but look for the final stationary configura- 
tions. Due to thermodynamical considerations we expect 
the condensate to achieve the least energy for given ex- 
perimental parameters, {N, g, f2}. Such configuration is 
called a ground state. 

Within the framework of our mean-field model, a sta- 
tionary state is represented by a wavefunction with the 
form 



(r). 



(4) 



The first way to find such solutions is to introduce (g) 
in our mean-field model ([!]) and directly solve the result- 
ing partial differential equation. This method has several 
disadvantages: (i) It is difficult due to the many degrees 
of freedom that it involves, (ii) this new equation is sat- 
isfied not only by the ground state but also by excited 
states, and (iii) there is no mathematical guarantee yet 
for the existence of solutions of this problem. 

The second way to characterize a stationary is solution 
is by studying the energy functional and use the fact that 
each stationary configuration is a critical point of the 
energy (g) for given {N,g,Ct}, 



dE 

dip 



W = o. 



(5) 



N 



Since we are indeed looking for the ground state, one 
should use a minimization procedure able to find the min- 
ima of the energy functional for given parameters. 
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C. Reshaping the energy functional 

As we stated above, our objective is to find the solu- 
tions "0^ which are the minima of the energy subject to 
the restriction J \ip\ 2 = N, for a fixed angular speed of the 
trap and for a given interaction. The existence of at least 
one minima or ground state for this variational problem 
has been proved elsewhere |Q, they are guaranteed to 
be stable and represent the energetically most favorable 
configurations for given N and £1. 

Such a problem is called a constrained optimization 
and we will refer to the solution, tp^, as a constrained 
minimum. The existence of constraints in an optimiza- 
tion problem poses serious difficulties to traditional de- 
scent methods, since it is difficult to design a efficient 
minimization algorithm which takes care of the con- 
straints at each step 14|. To take into account this re- 



striction one should use Lagrange's multipliers, i.e. just 
adding a fraction of the constraint w(/z, N) (also called a 
"penalizer" ) to the original functional 



F[%j}] =E[tI>\+w{h,N). 



(6) 



It is not difficult to show that F[ijj] and E[ip] have the 
same stationary states, and that any absolute or relative 
minimum of F[tp] is also a constrained minimum of E [ip). 

One expects that a wise choice of the penalizer will 
establish a one-to-one correspondence between the values 
of the chemical potential, /i, the absolute minima of F[tp] 
and the ground states of our condensate. The advantage 
of F [ip] over E[i/j] is that in the new functional there is 
no need to care for constraints: the value of N for the 
ground state is determined by the chemical potential, fi, 
which is fixed throughout the minimization process. 

The traditional choice for a Lagrange multiplier leads 
to the definition of the free energy 



T[iJj] =E[xj)} ~ fiN[ijj]. 



(7) 



However, this functional is not bounded below and thus 
ipp is at most a local minimum of T . Consequently, this 
functional cannot be used to characterize the state of the 
condensate. This is an interesting result, since one could 
be tempted to think that the traditional definitions of the 
thermodynamical potentials are suitable to characterize 
all physical systems. We must remember, however, that 
our condensate is being described by a mean-field model, 
a simplification of a more complex model, and as such 
what works for the full problem needs not work for the 
simpler one. 

In our search for suitable functionals we have found a 
simple one which we call the nonlinear free energy 



F[tl>] = E[tl>] + \ (N[ip] - A) 2 . 



(8) 



First and most important, it can be proven that F[tp] has 
at least one finite norm absolute minimum for each value 
of A [O] , and that each of those minima corresponds to 



a constrained minima of the energy, E[tp]. And second, 
there exists a simple and invertible relation between the 
usual thermodynamic variables — the chemical potential 
and the number of particles, (/i, N) — , and our new vari- 
ables — Lagrange's constant and the number of particles, 
(A,Af)- 



,i = N[i/> M ] - A. 



(9) 



An important feature of our new functional (g) is that 
it highly nonlinear with respect to ip. This poses no 
new difficulty, since our original equations (|l|) were al- 
ready nonlinear. Indeed Eq. (|) has proven to be the 
most natural choice for many other problems, such as the 
propagation of incoherently coupled laser beams through 
saturable media (lal. 



D. Optimal methods for minimization 

To minimize the nonlinear free energy (||) we follow 
Rcf . |l(| . First we choose the right function space, which 
in our case is the Sobolev space H 1 (R 3 ) = {ijj/ijj,'Vip £ 
L 2 } of functions which admit at least one spatial deriva- 
tive. This space is equipped with a scalar product 



(r)0(r) + W(r) ■ V0(r)] d n r. (10) 



and a norm ||^'| = {tp,ip). To obtain an explicit expres- 
sion for the gradient of the functional in H 1 we perform 
a first order expansion of F[ijj] around a trial state ip 



F[if> + ed] = F[if>] -I 

+ 0(£ 2 ). 
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(11) 



We have to turn this expression into something which 
looks like Frechet's definition of a derivative. This means 
that we have to find some d> such that 



Re 



S— X7S dE 



= Re 



WSV(j>] . (12) 



If we integrate by parts and impose that this equality 
be satisfied for all perturbations, 5, the problem has a 
formal solution which is given by a Lagrange equation 



(I - A) 6= — = - V — . 



(13) 



In consequence, our formal expression for the Sobolev 
gradient of F\^\ finally reads V S F = (1 - A^Vf, 
where VsF stands for the Sobolev gradient, VF is the 
ordinary one, and (1 — A) -1 represents the inverse of a 
linear and strictly positive definite operator. 

We may now use the Sobolev gradient of our functional 
as the direction of descent for a minimization procedure 
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-^A + y + g\u\ 2 -flL 2 



v. (14) 



The preceding equation converges to some stationary 
state </>(r) = lim T _ >00 ^(r, r). To grant convergence to 
the true ground state one needs to improve the descent 
method, using for instance the nonlinear conjugate gra- 
dient method instead of ( |l4| ) , and introducing some type 
of relaxation which helps avoid saddle node points. In 
any case, the minimization process must be performed 
on suitable space of functions that we built using a dis- 
crete Fourier basis with 64 x 64 x 128 modes. 
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FIG. 2. Three-dimensional surface plots of one to four 
vortex lines for e = 0.03 and Q/lj ± = 0.41,0.55,0.6,0.65. 
The walls of the boxes reflect different integrated views of 
the gas cloud (i.e. as they would look if imaged with high 
resolution), and within each box we show the lines of lowest 
density within the condensate. Each plot is represents a cen- 
tered box (x,y,z) 6 [—6,6] x [—6,6] x [—60,60], with sizes 
in adimensional units. Parameters are close to experimental 
values: Ng = 9000, e = 0.03, lu ± /uj z = 219/11.7. 



the number of atoms was chosen to match, Ng — 9000, 
which corresponds to a few times 10 5 Rb bosons but our 
results remain qualitatively valid for an ample range of 
gN values as will be shown below. 

The geometry of the trap is a very elongated one, 
u>± > 18.7lj z , and for the small transverse deformation 
of the trap we have tried e = 0.0 as well as e = 0.03, the 
latter being the closest one to the actual experiment ||. 
This means that we have studied both axially symmet- 
ric setups and completely asymmetric setups, and in both 
cases our study has provided essentially the same results. 
This is due to the fact that an intense nonlinear interac- 
tion (Ng ~ 10 4 ) effectively prolongs the existence of all 
solutions from the symmetric setup to the asymmetric 
setup fPj\ . 

The first result is that an increase of the angular speed 
causes the minimum of the energy to move from a node- 
less ground state into states with one, two, three and more 
vortices. Consequently, as it was already predicted 0, 
there exists a cascade of increasing angular frequencies, 
{fli < SI2 < . . .} for the nucleation of one, two and more 
vortices — the larger the rotation frequency, the more 
vortices, and for our setup, fli ~ 0.4(0), — 0.5(3), 
fl 3 ~ 0.5(8), and {I4 ~ 0.6(4). Figure |l| shows the struc- 
ture of a ground state hosting two vortices at £1 = 0.55. 

The second and most important result, which is al- 
ready evident in Fig. [l] is that vortices are nucleated 
with a stable and stationary bent shape even in axially 
symmetric traps. In Fig. we show three-dimensional 
pictures of a condensate with up to four bent vortex lines. 

These contorted shapes are absolutely stable configura- 
tions which lack rotational symmetry of any kind, even 
discrete, something which was thought to be forbidden 
according to Ref . . In that case asymmetric states were 
found only for discrete values of the angular speed, right 
on the transition between different number of vortices, 
and always exhibiting energetic instabilities. Further- 
more, due to the bending, while the trunk of the conden- 
sate rotates, the caps remain almost still. This feature 
allows the gas to accommodate a fractional value of the 
angular momentum which is between and 1, and which 
never reaches 1, as it was already conjectured in [|l(|. As 
the speed of the trap increases, the angular momentum 
grows continuously by means of pulling the vortex line 
straighter and not only with discontinuous jumps. 



III. RESULTS 

A. Bent vortex lattices in current experiments 

We have applied the numerical methods outlined above 
to different setups. The family of numerical experiments 
in this subsection resembles the experiments of Madison 
et al. ||||. For these experiments with 87 Rb conden- 
sates, the bosonic interaction is ruled by the scattering 
length as — 5.5 nm. For several pictures in this paper, 



B. Regularity of the vortex lattice 

From old studies with liquid helium and more recent 
works with condensates of simpler geometry, it has been 
expected that vortices should form regular triangular lat- 
tices. On the other hand, our finding above seems to 
suggest that the vortex aggregate rather adopts an ir- 
regular shape. As we will show below, the condensate 
actually develops a medium range order, which dilutes in 
the vicinity of the cloud boundary. 
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In Fig. |3| we use the solutions of our previous section to 
simulate the pictures that should be seen in experiments. 
The first row of two-dimensional plots [Fig. ||(a-d)] are 
upper views of the condensate, and represent the col- 
umn density of the bosonic cloud along the transverse 
directions. These pictures bear an extreme resemblance 
with the experimental photographs of Refs. Hf|, show- 
ing blurred clouds where the vortex cores seem partially 
filled. 

In Fig. ||(i-l) we show the density of the conden- 
sate as seen at half the height of the condensate, i.e. 
\ip{x,y,z = 0)| 2 . At the moment of revising this paper 
an experiment has been made which observes this type 
of transverse cuts of a stirred sodium condensate . In 
that experiment a thin slice of bosons is promoted to a 
different internal state and then imaged separately from 
the rest of the cloud. The resulting pictures are like those 
in Fig. ||(i-l), where the regularity of the Abrikosov vor- 
tex lattice is made evident. The blurring has disappeared 
and the holes arrange along the expected triangular lat- 
tice. In the above mentioned paper the bending of the 
vortex lines is also reported. 

As further confirmation we suggest that using the 
setup of Ref. ||] and watching the condensate not from 
above, but from one side would allow a direct observa- 
tion of bent vortex lines lead to pictures as those in Fig. 
(e-h). 
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FIG. 3. Condensate shapes for the same states as in Fig. 
f§ Q/Vl = 0.41 (a,e,i) ,Q/uj ± = 0.55 (b,f,j), Q/w± = 0.6 
(c,g,k) and f2/u>j_ = 0.65 (d,h,l). Shownare up (a-d) and side 
views (e-h) and (i-1) two-dimensional cuts of the cloud at half 
its height. 



C. What is the ultimate cause of bending? 

The fact that the bending of the vortex line takes place 
not only on the asymmetric trap of current experiments, 



but also on an radially symmetric trap, e = represents 
a surprising novel type of symmetry breaking in which 
the superfluid not only chooses the sense of rotation, but 
also a plane for its bending. 

It is important to emphasize the counterintuitive na- 
ture of the vortex bending in the case of an isolated vor- 
tex line. First, opposite to the case of 4 He, the BEC is 
not constrained by any recipient and there are no asym- 
metric boundary conditions which could easily explain 
the deformation of the fundamental solutions. Second, 
although qualitative reasons for such bending may be 
found a-posteriori, they have never been reported be- 
fore. And even though physical arguments may justify 
the bending of the vortex line, they can hardly support 
the fact that the bent vortex line is a stable stationary 
configuration. Both the stability and the stationarity can 
only be demonstrated with a direct study of the energy 
functional. 
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FIG. 4. Side views of the condensate for N = 10 5 atoms 
of 87 Rb, and decreasing elongation. The trap parameters and 
plot dimensions (fi, u>±/2tt, uj z /2tt, r max , z max ) are (a) (0.5wx, 
219 Hz, 11.7 Hz, 2.63 nm, 49.09 nm), (b) (0Alu ± , 96 Hz, 25 
Hz, 5.00 nm, 19.23 nm), and (c) (0.3w_l, 50 Hz, 50 Hz, 8.51 
nm, 8.51 nm). (d) Uses the same trap as in (a), but with 
{iV = 10 4 , Q = 0.53, r max = 1.65, z max = 30.98}. 

Nevertheless, it is legitimate to ask for the ultimate 
reasons of this rupture of symmetry. Since our search of 
states was based on energetical considerations, we have 
studied the dependence of the bending with respect to 
each of the free parameters in our problem: N, Q, g, e, 
and the elongation 7 = lu±/lo z . 

First, the influence of e is discarded. As we mentioned 
above, bending exists in either symmetric or asymmetric 
traps. Rather, the asymmetry seems to take part only in 
the dynamics of the nucleation process, by changing the 
value of the critical frequencies and inducing a type of 
hysteresis [p"of . 

A different role is played by the elongation of the trap, 
7, and the effective interaction U = Ng. Starting with 
a configuration and lowering either the elongation or 
the number of atoms we see that the bending becomes 
smaller and eventually disappears [Fig. [|(a-c)]. 

We have tried to measure the bending so as to deter- 
mine the minimal enlogation and interaction which are 
required to induce this phenomenon. The problem is that 
due to the small changes of L z , the full minimization pro- 
cedure does not allow us to bound this minimal elonga- 
tion accurately. An alternative procedure is to study the 
relation between the angular speed which is required to 
stabilize a straight vortex, Qi and the angular speed at 
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which the ground state acquires some angular momen- 
tum, Qi. The first value arises from the study of normal 
modes around a straight vortex, while the second value 
is the energy difference between a straight vortex and 
a vortexless state [[To). When Cli ^> Qi there are val- 
ues of the angular speed where bending can be favorable. 
This phenomenon has been referred to in the literature as 
"anomalous modes" 1 8 , and it is a signature of bending. 

In Fig. H we plot the difference f2i — fli, as a func- 
tion of the elongation, 7, and of the effective interaction, 
U = Ng for an axially symmetric trap. According to this 
study, a large elongation of the trap is required to make 
a bent vortex energetically favorable over a straight one. 



tice at half the height of the condensate and deform not 
very far from the condensate core. 
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FIG. 5. Grayscale plot and contour lines for the difference 
Oi — fli between the speed required to nucleate a vortex, fii, 
and the speed required to stabilize a straight one, fii, as a 
function of the elongation of the trap, 7 = lu±/uj z , and the 
effective interaction, U = Ng. A horizontal dash-dot line 
marks the spherically symmetric trap, 7 = 1. For Q\ > Qi, 
i.e. above the thick solid line, the bending of vortices becomes 
energetically favorable. All figures are adimensional. 



IV. CONCLUSIONS 

In this work we have studied vortex lines in a very 
elongated Bose-Einstein condensate in a rotating trap. 
By d eveloping both a customized energy functio nal (S ect. 
1IC) and an optimized descent method (Sect. 



II D) we 



have been able to find the ground state of a condensate 
in a family of different experimenta l confi gurations. 

The first part of our study (Sect. Ill A ) focused on re- 
alistic values of the experimental parameters, taken from 
the work by Madison et al. [§JiJj]. The 

main conclusion 

of this work is that by increasing the rotation speed, the 
condensate achieves states with one, two, three and more 
bent vortices. Such vortices form a regular Abrikosov lat- 



Appart from the observation reported in 12 and while 
current experiments at the ENS [|,||,[l9| are carried on a 
regime in which vortex lines should bend, this has not 
been yet observed there. The reason is that the imag- 
ing of a condensate which is a few fj,m large requires 
a previous expansion in which the cloud becomes disk 
shaped and bending cannot be appreciated. Further- 
more, the study of the elongated cloud without vortices 



|10,|20| shows that the nucleation of vortices is subject 
to hysteresis and one must actually exceed some rotation 
frequency fl m , which is rather large and makes it difficult 
to selectively produce one, two or more vortices. 

On the other hand the bending of vortex lines has sig- 
natures which are also observed in current experiments. 
First, the bending by itself can explain the apparent fill- 
ing of vortex lines when seen from above [EJ , and the reg- 
ular pictures which appear in most recent observations 
with sliced sodium condensates |Q . Second the bending 
also accounts for the continuous growth of the angular 
momentum with respect to the rotation speed, and the 
fractional values of the angular momentum, < I < 1, 
which arise in the indirect measures of [§ and g9|. And 
finally, the nucleation of many bent vortices and their 
subsequent interaction during the expansion phase may 
lead to turbulent structures that should explain the lack 
of regular pictures above a certain angular speed J|,|| ■ 

The second part of this paper studied the dependence 
of the bending with respect to the parameters of the mean 
field model ([!]). Here we conclude that it is both the 
elongation and interactions which induce the bending of 
the vortex line, while the transverse asymmetry plays no 
important role. 

Our results imply that past studies devoted to the 
quasilinear limit (U <C 1) would become of no applica- 
bility for elongated traps due to the lack of bending in 
the simplified models. Second, along this line it would be 
nice to find analytically the geometry of the vortex line 
as a function of the interaction U and the elongation 7. 

Our study also reveals that current experiments are 
being developed in a regime which is qualitatively differ- 
ent from the setups that have been studied up to date 
Values of the relevant parameters (7, U) in cur- 
rent experiments are so far from most studies that new 
methods must be developed to accurately describe these 
amazing systems. 

As an example, let us pose one of the technically diffi- 
cult questions that arise in this work. From our numeri- 
cal work it seems that the extremes of some of the vortex 
lines actually reconnect not too far away from the core of 
the cloud, forming what is called a vortex ring. Although 
it is not possible with current numerical methods to fully 
support this conjecture, vortex rings have already been 
observed in BEC . The limits of zero and infinite ra- 
dius of a vortex ring lead to a dimcnsionless zero and an 
isolated vortex line, respectively. Therefore vortex rings 
would be a nice tool for explaining the nucleation of vor- 
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tices and they would allow to interpret the states in this 
work as the result of an incomplete nucleation. The proof 
or refutation of this conjecture remains a mathematically 
challenging problem. 

This work has been supported by grant BFM2000- 
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